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Abstract 

^sO | We show that, for any n > 0, the Heisenberg interaction among 2n qubits (as spin-1/2 

particles) can be used to exactly implement an n-qubit parity gate, which is equivalent in 
constant depth to an n-qubit fanout gate. Either isotropic or nonisotropic versions of the 
interaction can be used. We generalize our basic results by showing that any Hamiltonian 
\ (acting on suitably encoded logical qubits), whose eigenvalues depend quadratically on the 

Hamming weight of the logical qubit values, can be used to implement generalized Mod 9 gates 
for any q > 2. 



This paper is a sequel to quant-ph/0309163, and resolves a question left open in that paper. 

1 Introduction 



O 
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Let TC be the Hilbert space of n qubits, where n > 1. The fanout operator F n : H —* H, depicted 
in Figure ^ copies the (classical) value of a single qubit to n other qubits. Unbounded fanout is 
usually taken for granted in models of classical Boolean circuits, even when the fanin of gates is 
bounded. They cannot be taken for granted in quantum circuits, however, since copying the value of 
a quantum bit to n— 1 other bits requires significant nonlocal interactions. Unbounded fanout gates 
have been shown to be a surprisingly powerful primitive for quantum computation, allowing one to 
reduce the depth of a circuit computing, say, the Quantum Fourier Transform (QFT) to essentially 
constant depth JHJ- The quantum part of Shor's factoring algorithm can thus be implemented in 
constant depth if unbounded fanout gates are available. This result about the power of fanout is 
especially important because in most of the significant proposals for implementing quantum circuits 
to date, long computations will surely be difficult to maintain due to decoherence, current quantum 
error correction techniques notwithstanding. Shallow quantum circuits may prove to be, at least in 
the short term, the only realistic model of feasible quantum computation, and fanout gates would 
increase their power significantly. 

Without some quantum gate with unbounded width (arity), it is not clear that any nontrivial 
decision problem can be computed by o(log n)-depth quantum circuits with bounded error. This is 
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Figure 1: Definition of the fanout gate. 



certainly true if we only allow a one-qubit output measurement (then the output can only depend 
on 2°( logn ) input qubits — see [3] for a discussion), but it also seems to be true even if we allow a 
computational-basis measurement of an arbitrary number of qubits at the end Even if we allow 
unbounded AND gates (generalized Toffoli gates), it is not clear what we can do in sublogarithmic 
depth. We do know that we cannot approximate fanout gates this way jSJ. 

To summarize, fanout gates are an extremely useful, and perhaps necessary, primitive for al- 
lowing small-depth quantum circuits to solve useful problems. Furthermore, implementing fanout 
with a conventional quantum circuit requires logarithmic depth, even if unbounded AND gates are 
allowed. Therefore, implementing a fast fanout gate will require an unconventional approach. 

We provide such an approach here by showing that the fanout operator arises easily by evolving 
qubits via a simple and well-studied Hamiltonian, the spin-exchange or Heisenberg interaction, 
together with a modest amount of encoding and decoding of qubits (which only requires constant- 
width gates and constant depth). Our results answer positively a question by I. L. Chuang, who 
asked how certain forms of the Heisenberg interaction, which are implementable in the laboratory, 
may be useful for quantum computation [31 H]- In particular, we show that the fanout gate on n 
logical qubits can be achieved exactly by encoding them into 2n physical qubits (each a spin-1/2 
particle), then applying the Heisenberg interaction to the encoded qubits. The interaction need 
not be isotropic; both isotropic and nonisotropic versions of the interaction work equally well. 

In [7j, we showed that a variant of the Heisenberg interaction, where the Hamiltonian is pro- 
portional to the square of the z-component of the total spin, can implement parity easily (without 
encoding). When applied to three qubits, this interaction yields an "inversion on equality gate" 
/=, defined by 

/= i \ = J ~\ X V Z ) tfx = y = z, 
\ \xyz) otherwise. 

1= and single-qubit gates together form a universal set of gates. Recently, implementation of 1= 
as well as the three-qubit parity and fanout gates in NMR using the above Hamiltonian has been 
reported (SJ. Our current paper affirmatively answers a question left open in [2] as to whether par- 
ity/fanout can be implemented using more common forms of the Heisenberg interaction, involving 
x-, y-, and z-components of the total spin. 
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Figure 2: Definition of the parity gate. 



In Section [2] we define the general Heisenberg interaction between n identical spins, as well 
as the special case of interest to us. In Section 13.11 we give an implementation of the (r + l)-bit 
parity gate, depicted in Figure EJ where r = n/2. The fanout gate arises by placing Hadamard 
gates on each qubit on both sides of the parity gate (2(r + 1) Hadamard gates in all), and thus 
implementing parity is equivalent to implementing fanout. We also show in Section l3.2l how different 
qubit encoding schemes can reduce the ratio n/r to be arbitrarily close to one. In Section 0J we 
generalize our results in two ways: (1) any Hamiltonian whose eigenvalues depend quadratically on 
the Hamming weight of the logical qubits can be used to implement parity, and hence fanout, and 
(2) any such Hamiltonian can implement generalized Modq-gates directly for any q > 2. 



2 Preliminaries 

The Heisenberg interaction describes the way particles in the same general location affect each 
other by the magnetic moments arising from their spin angular momenta. Given a system of m 
identical labeled spins described by vector operators Si, ... , S m , the Hamiltonian E of the system 
is a weighted sum of the energies of all the pairwise interactions, plus a term for any external 
magnetic field (assumed to be in the z-direction): 

E = -^2Ji, j S i -S j + aY J (S l )z, (1) 

i<j i 

where the Jij and a are constants. 1 In this paper, we will show how this interaction can implement 
fanout in the special case where all the Jij are equal. In this case, E is related to the squared 
magnitude of the total spin of the system. 

We will assume here that physical qubits are implemented as spin-1/2 particles, with |0) being 
the spin-up state (in the positive z-direction) and |1) being the spin-down state (in the negative 

1 The Jij are usually assumed to be positive, appropriate for ferromagnetic interactions which give the lowest 
energy when spins are aligned in parallel. The value of a is the product of the magnetic field strength and the 
gyromagnetic ratio for the individual spins (see |l 21 §21.3] for example). 
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z-direction). Given a system of m qubits labeled 1, . . . , m, we define 



1=1 



i=l 




i=l 



where X{, Yi, and Zi are the three Pauli operators acting on the i'th qubit. J x , J y , and J z give 
the total spin in the x-, y-, and z-directions, respectively. The squared magnitude of the total spin 
angular momentum of the system is given by the observable 

j j j. | J y | J g • 

Note that 

j2 = — I+ 2 E (W J+ Y l Y J + Z l Z J ) = —I + Y,Si-S 3 , 

l<i<j<m i<j 

where Si = \{Xi, Yi, Zi) is the vector observable giving the spin of the ith qubit. It is then clear 
that, in the absence of an external magnetic field, J 2 is linearly related to the energy E above. 
This is an isotropic Heisenberg interaction. To account for an external field in the z-direction, we 
define, for any real a, 

H a = -J 2 + aJ z . (2) 

This is the case of the Heisenberg interaction where all the Jj are unity. 2 Our methods can also 
accommodate an extra term in the Hamiltonian proportional to J 2 with no additional effort (see 
Section [3. so we consider the more general Hamiltonian 

Ha,/3 = -J 2 + aJ z + f3J 2 z (3) 

for any real a and (5 such that (3 ^ 1. We will evolve the system of m qubits using H a ^ as the 
Hamiltonian. A related Hamiltonian J 2 is used in [7] to implement fanout; this is an easy case, 
since each computational basis state is already an eigenstate of J z and thus the implementation 
requires no encoding of qubits. Using the current Hamiltonian H a ^ is more complicated and 
requires encoding logical qubits into groups of physical qubits, so that the tensor product of all the 
physical qubits encoding a logical basis state will be an eigenstate of H a p. 

In the sequel, we choose units so that H = 1. If A and B are both vectors or both operators, we 
say, "A ex B" to mean that A = e B for some real 9, that is, A = B up to an overall phase factor. 
We use the same notation with individual components of A and B, meaning that the phase factor 
is independent of which component we choose. If x £ {0, 1}™ is a bit vector, we let wt(x) denote 
the Hamming weight of x, that is, the number of Is in x. 

2 This does not lose generality, since constant factors in the energy can be absorbed by adjusting the time of the 
interaction. 
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2.1 Spin States 

The properties of the J-operators are well-known. See, for example, Bohm 2 . We will review 
the essential ones here. The commutation relations are [J x , J y ] = i,J z , and likewise for the two 
other cyclic shifts of the indices. J 2 commutes with J z , so one may choose an orthonormal basis 
of the n-qubit Hilbert space TC that diagonalizes both simultaneously. Eigenstates of J 2 and J z 
are traditionally labeled as \j,m,£), where < j < n/2 and —j<m<j, with n/2 — j and 
j — m both integers. We have J 2 |j, m, £) = m\j, m, £), and J 2 \j, m, £) = j(j + 1)| j, m, £). The extra 
parameter £ is used to give distinct labels to different basis vectors in degenerate eigenspaces of J 2 
and J z . These basis vectors can be chosen so that, for any value £ that appears as the third label of 
some basis vector, the basis vectors labeled by £ span an irreducible spin representation, that is, a 
minimal subspace of TC invariant under the action of J x , J y , and J z . This space will be spanned by 
the basis vectors \j, -j,£), \j, -j + !,£),..., - l,£),\j,j,£), for some j = j(£) depending only 
on the label £, and is called a spin-j representation. Letting J + = J x + iJ y and J_ = = J x — i,J y 
be the usual raising and lowering operators, respectively, we may adjust the phases of the basis 
vectors so that 

J+\j,rn,£) = + 1) -m(m+ 1) \j,m + 1,1), 

J-\j,m,£) = + 1) - m(m - 1) \j,m- 1,£). 

(This sets the relative phases of states within the representation, but still allows the overall phase 
of the representation to be adjusted relative to other representations.) 

2.2 Number of Spin Representations 

An important fact that we will use later is that for each j, there are exactly k n j := ( n /2-j) — 
(n/2-j-i) man y spin-j representations in the decomposition of Ji? One way to see this is as follows. 
For any j > such that n/2 — j is an integer, let TCj be the eigenspace of J z with eigenvalue j (if 
j > n/2, then Tij has dimension zero). Clearly, dim(Tij) = \ n i2-j)i since TLj is spanned by all the 
computational basis vectors with Hamming weight n/2 — j. The J + operator maps TCj into Tij+i, 
and so its kernel on Tij has dimension at least k n j > 0, given above. Now the space ker(J + ) n TLj is 
spanned by the set of all states of the form \ j,j, £}, so there are no less than k n j distinct values for 
£ occuring in the set, and each one labels a distinct spin-j representation. Finally, since each spin-j 
representation has dimension 2j + 1 for all j, and since TC has 2 n dimensions, a simple counting 
argument shows that there can be no more than k n j many spin-j representations, either. 

It follows that there are a total of (^ n " 2 j) many spin representations in the decomposition of 
TC. Moreover, if n is even, then it is easy to show that the representations are evenly split between 
those where n/2 — j is even and those where n/2 — j is odd: ^( n " 2 ) representations for each. This 
fact will be used in Section [3.21 

2.3 Spin States versus Computational Basis States 

Finally, we mention how some of the spin states relate to computational basis states. There is one 
spin-n/2 representation in the decomposition, namely, the completely symmetric representation, 

3 By convention, if k < then (™) = 0. 
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Figure 3: A two-qubit encoder. 



which is spanned by the states 

\n/2,n/2-k)=(fj ^ £ \x) 

^ ' wt(x)=k 

for integer k with < k < n. (This equation sets the overall phase of the spin-n/2 representation.) 
A key point in this paper is to note that |§, §) = |0 n ) is a tensor product of single qubits. This 
means that some spin states involve little or no entanglement among qubits and thus can be prepared 
using only reasonably local interactions. More generally, suppose we group some of the n qubits into 
disjoint pairs (it, ji), (i 2 , 32), (ip,3 P ) for some < p < n/2 (that is, . . .,i p ,j p G {1, . . . , n} 

and are all pairwise distinct), then we form a state by putting each pair (ik,jk) of qubits into 
the singlet state 1 0, 0) = (|10) — |01))/\/2 and each of the rest of the (unpaired) qubits into the |0) 

S t ctt 6 1 6 

where j = n/2 — p, and S is the set of all unpaired qubits. Then it is easy to check that \ip) is an 
eigenstate of J z (with eigenvalue n/2 — p) and is in ker(J + ). It follows that is also an eigenstate 
of J 2 , because = J-J + \ip) = (J 2 — J 2 — J z )\tp). Many of the states \j,j,£) where j = n/2 —p can 
be defined this way, but not all — different choices of the p pairs do not always produce states that 
are orthogonal to each other, or even linearly independent. 



3 Main Results 

3.1 Parity Gate by Heisenberg Interactions 

Recall the Hamiltonian H a ^ of © on the space TL of n qubits labeled 1, . . . ,n. H a ^ commutes 
with J 2 , so it has eigenvectors |j, m,£) with respective eigenvalues —3(3 + 1) + am + (5m 2 . 

Let x = x\ ■ ■ ■ x r be a vector of r bits, where r is some number no greater than n. We wish 
to encode the r-qubit computational basis state \x) into an n-qubit eigenstate \j, m,£) of so 
that 3 depends linearly on wt(x), and we want to do this by using gates that act on as few qubits 
as possible. The easiest way to accomplish this is to have n = 2r and create encoded states of the 
form given by ©. One may encode each input qubit (with an ancilla) into two qubits, sending 1 00) 
to |0l) := |00) and sending 1 10) to \lz) ■= ( 1 10) — |01))/\/2, the singlet state. A simple circuit — one 
of many — for this is shown in Figure 0J which defines the encoding operator E. (Since = E, we 
will also use E to decode.) Clearly, there are other operators that will work just as well, since E is 
underdetermined. We encode the ith input with ancilla into the qubits 2i — 1 and 2i oiTL. Thus 
\x\ ■ ■ ■ x r ) maps to \xl) := <8> • • • &> \x r ,L), and this state is in the form of (HJ), where the set S 

consists of all physical qubits encoding the |0l) states, i.e., S = Uz i= o{^ — 1 , 2z} . If x ^ y, then 
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clearly (xl\Vl) = 0, so we can assume without loss of generality that \x£) = \j(x),j(x),x), where 
j(x) := n/2 — wt(x), and x itself is used for the label. 
Suppose that (3^=1, and let 



t :-- 



7T 



2|/?-l| 



We let 



U : = e - itHa -e, 



(5) 



(6) 



the unitary operator resulting from evolving the qubits with H a ^ for time t. For fixed input vector 
x = x\ ■ ■ ■ x r € {0, l} r , let k = wt(x), and for b G {0, 1} let x b := x\ ■ ■ ■ x r _\b. To compute the 
parity of k with a quantum circuit on input \x), we first run qubit r through a Hadamard gate to 
produce the state 

\x°) + (-l) Xr \x l ) 



\<Px) ■-- 



V2 



Next, we encode each qubit as described above to obtain 

|x° ) + (-!)*'• |4) |j ,jo,x°) + (-1)^^, ji.x 1 ) 

M := 71 = 71 ' 

where we have set % := j(x b ) = n/2 — wt(x fc ). Observing that j\ = jo — 1, we see that \*p x ,L) is 
a balanced superposition of two spin states, one with even j and the other with odd j. We now 
apply U to \<p x , L )- Noting that H a ^\j,j,£) = - l)j 2 + (a - l)j)|j, j,^) for any j, we set 

a — 1 



7 := 



0-1' 



and define 



l^> := U\j,j,£) = e- is ^ + ^ 2 \j,j,e) 
for any j and £, where for convenience, we are letting 



(7) 
(8) 



s :-- 



1/3-11 



1 if /? > 1, 
-1 if/3<l. 



We have 



= (hio,xO> + (-l) Xr |^ ll xi»/V2 

= [exp(-isvr(jo + 7Jo)/2)|jo, jo,^ } 

+ (-l) a! "exp(-i a 7r(j 1 2 +7ii)/2)|ji,ji,x 1 >] /72 
= |xi,l) (8) ■ ■ ■ <8> |x r _i )L ) (g) |* x ), 



where 



,-isir(j2+yj )/2 



2 |0 L ) + (-l)^ e -«^i+^)/ 2 |l L ) 



72 



-«s7r(jg+7j )/2 



72 



oc 



1 

71 
i 

TIL 



|o L ) + (-i)^ e - iCT (- 2j0+1 - 7 )/ 2 |i L ) 

|OL) + (-l) io+a; '-e^ (7 - 1)/2 |lL) 
|0 L ) + (-lJ^We^-^lli) 



(9) 

(10) 

(11) 
(12) 
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Figure 4: Circuit to implement parity with Heisenberg interactions. 



is the state of the (2r — l)st and (2r)th qubits. (Recall that j\ = jo — 1, and that jo = § — wt(x°) = 
r — wt(x) + x r .) 

For any y = y\ ■ ■ ■ y r £ {0, l} r , it is easy to see that if wt(x) and wt(y) have opposite parity, 
then \^ y ) and |^ x ) are orthogonal: it follows immediately from (|^V ll2() that 

2(* w |* x > oc 1 + (-l) wt (*)+ wt (») J 

which is zero if wt(x) and wt(y) have opposite parity. This analysis shows that we have isolated the 
parity information in the rth logical qubit. We decode just the two physical qubits corresponding 
to this qubit to obtain the state 



-is7r(jg+7jo)/2 r 



V2 



|00) + (-l) r + wt W e ^(T-l)/2|l ) 



1 

V2 



|0) + (-l) r + wt ( a; ) e is7r (7-l)/2| 1 )" 



The second of these qubits is the restored ancilla. We then apply two gates, V followed by H, to 
the first qubit, where 

" 1 

e -is7r(2r+7-l)/2 



V 



and H is the Hadamard transform. This yields |wt(x) mod 2) as the state of the first qubit, up to 
some unconditional phase factor. 

If we do not mind the extra phase factor, which only depends on 7, r, s, and wt(xo), then 
we can simply use E to decode the other pairs of physical qubits, and we thus obtain a circuit 
that computes parity. To cleanly and exactly match the parity gate defined in Figure El however, 
we may first copy the parity information onto a fresh qubit, then undo the previous computation. 
The circuit for the latter operation is shown in Figure 0] The gate on the left in Figure E] is an 
(r + l)-qubit parity gate. 
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Figure 5: A gate that encodes c logical qubits x = x\ - ■ ■ x c into d physical qubits in state \j x ,jx,£x)- 



Remark. The circuit of Figure 0] uses the gate. Any unitary gate U' that agrees with IP on 
the subspace TL' of H spanned by vectors of the form \j,j,£) can substitute for U> to implement 
parity exactly. We would like to implement U' by evolving the same Hamiltonian H a g for some 
positive length of time. We can do this if there is a u > such that e~ luHa ' 13 fixes all vectors in 
H' . Then we can implement U' by evolving H a a for time ku — t where t is given by © and k 
is some integer such that ku > t, i.e., U' := e~ t ^ ku ~ t ^ Ha ' f3 . By 0, it can be shown that such a 
u exists if and only if the 7 of (0) is rational. For arbitrary real 7, we can still implement U' by 
evolving the altered Hamiltonian H a i pi for some suitable ot',/3 1 and time tf. Assuming /3' 7^ 1 and 
letting 7' := %rE% 1 h follows from (JSJ) by a straightforward argument that for t' > 0, the operator 

U' := e~ ltHa ''i 3 ' is a suitable replacement for W (i.e., f7'C/ fixes all vectors in 7i') if and only if 
there is an integer i such that (i) (2£ + 1)(7' + 1) + s(7 + 1) is an integer multiple of four, (ii) 21 + 1 
and — 1 have the same sign, and (hi) t' = 2(^-1) • 

3.2 More Compressed Encodings 

In the previous section we encoded each logical qubit into two physical qubits before applying U. 
By encoding groups of logical qubits, we can reduce the physical-to-logical qubit number ratio as 
close as we want to unity. Fix an integer c > 0, and let d be the smallest even integer such that 
{J/2) — 2 C - We can compute parity as before by dividing the logical qubits into groups of c qubits 
each (assume for convenience that r is a multiple of c), and encoding each group into a group of d 
physical qubits, yielding a ratio of d/c. Let E' be such an encoder, depicted in Figure 00 Our only 
requirement for E' is that it map each state |xO d-c ), where x G {0, l} 6 , into a state of the form 
\jx,jx,£x) where j x has the same parity as d/2 — wt(x) and l x ^ £ y if x ^ y. By our discussion in 
Section [2.21 there are enough spin representations on d qubits to allow this, so such an E' exists. 

Now by the considerations of Section \2. 31 we see that each input basis state \x) with x 6 {0, l} r 
is thus encoded into a superposition of spin states ai\jx, ji,£x) + a 2 \j2, j2, £2) + • • • , where all the ji 
are integers with parity equal to that of n/2 — wt(x). We then see by linearity that we can simulate 
the parity gate with a circuit identical to that shown in Figure HJ except that E is replaced with 
E' or E'^ as appropriate, and each encoding group has d physical qubits. 
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Since ( d / 2 ) = 2 d / y/d, we see that 

d _ d ^ 
c d—\ log 2 d 

which shows the trade-off between the size of E' and the ratio d/c. 

Encoding with an odd number of physical qubits per group is also possible, and may sometimes 
lead to a slightly better trade-off. For example, there are enough spin representations on five 
physical qubits to encode three logical qubits. 



4 Generalized Mod Gates from Any Quadratic Hamiltonian 

In this section, we show how to implement a Mod g gate directly, for any q > 2, using any Hamilto- 
nian whose eigenvalues depend quadratically on the Hamming weights of the inputs. More specifi- 
cally, we assume a Hamiltonian G n acting on n qubits, real constants a n , b n ,c n with a n > 0, and an 
encoding procedure E such that, for any computational basis state \x) over an appropriate number 
of qubits, S(|x)|00 • • • 0)) is an eigenstate of G n with eigenvalue a n w 2 + b n w + c n , where w = wt(x) 
and |00 - - - 0) is some ancilla state. Under these assumptions, we construct circuits implementing 
Modq gates for any constant q > 2, using evolution under G n . It is already known that the Mod g 
gates for all q > 2 are constant-depth equivalent to each other [Sj, so in effect, we already have a 
Hamiltonian simulation of any Mod g gate, via our implementation of the parity (Mod2) gate and 
the simulation in 9 . Our approach here is much more direct, however. Furthermore, it is only 
marginally more difficult conceptually to generalize our simulation to all q, rather than just q = 2. 
Our present development also subsumes the results in [7], where we implemented parity {q = 2) 
using the Hamiltonian J 2 . 

Fix q > 2. We consider q to be constant. The Mod g gate is a classical gate that acts on r 
control bits and a target bit. The target bit is flipped iff the Hamming weight of the control bits is 
not a multiple of q. We will actually simulate a more powerful version of this gate, the generalized 
Mod 9 gate, which has r control bits and q — 1 target bits £]_,..., t q — \. If w is the Hamming weight 
of the control bits, then the target bits t\, . . . , are all flipped, where i = w mod q, and the other 
target bits are left alone. Figure shows how to simulate a (standard) Mod 9 gate with a circuit 
using two generalized Mod g gates and a CNOT gate. 

We use some G n to implement the generalized Mod q gate via the circuit shown in Figure [7| 
Given an initial basis state \x) = \x\-- ■ x r ) of the control qubits, we first prepare q — 1 ancilla 
qubits into a state 

9-1 
3=0 

where the Cj are any fixed scalars such that \cj\ = 1/^/q (we may take Cj = 1/^fq for all j, for 
example). Note that |x)|y?) is a superposition of basis states with respective Hamming weights 
wt(x), wt(x) + 1, . . . , wt(x) + q — 1. By assumption, we have an encoder E that maps each compu- 
tational basis state \y) with y £ {0, 1} T+9_1 (possibly with additional ancillas) to a state \yi) over 
some number n of qubits such that 

c n) \Vl) ■ 
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Figure 6: Simulating a standard Mod g gate using generalized Mod g gates. There are r control 
qubits, and the ancillas on the right are the qubits labeled t±,t2, ■ ■ ■ ,t q -i, i.e., the target qubits of 
the generalized Mod g gates. 




W) 



Figure 7: Implementing a generalized Mod g gate using the Hamiltonian G n . Here, U = e~ ltGn for 
an appropriate t > 0, and R is described below. Any extra ancillag used by the encoder E are not 
shown. 
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Thus E maps to the state 



9-1 

3=0 



Next we apply U := e ltGn to \i[>x,l), where t = and k > is some fixed integer that is prime 



qa n 



to q (we may take k = 1, for example). Letting u; := wt(x) and 6 := b n /a n and c := c n /a n , we have 

9-1 

l/|^ SiL ) =^ Ci e- i7rfc [^') 2 + 6 ^+^+ c ]/9|(a;l J '0 9 - 1 - J ') L >. 
j=o 

We decode this state using E* to obtain the state |x)| 1 I / , J ,), where 

9-1 

|vj/^ ._ ^ c . e -iKk[(w+j) 2 +b(w+j)+c]/q\^3Qq-~t-j^ 
3=0 

is the state of the q—1 ancilla qubits. 

To see that we have isolated the value w mod q in the ancillae, we need only check that, for 
any integer v, (^ v \^ w ) = if v ^ w (mod q). Note that all the states of the form \^ u ) lie in 
a (/-dimensional subspace TL" of H, spanned by { FO 9-1- ^) | < j < q\. Let v be the Hamming 
weight of z. Then we have 

9-1 

(*«|* tt > = ^ l c i! 2 exp ( i7r M( v + J') 2 + + i) + c - (w + j) 2 - b(w + j) - c]/q) 
3=0 

9-1 

= q^ 1 ^2 ex P (^7r^[w 2 — w 2 + 2j(^ — + 6(u — w)]/q) 
3=1 
9-1 

cx q~ l exp (2iirjk(v — w)/q) 

3=1 

0(v mod q),(w mod g); 

where <5 Xj y is the Kronecker delta. Thus there is an orthonormal basis {\ctj) | < j < q} for 7i" 
such that, for all integers w > 0, there are real values 9 W such that \^ w ) = e l9w \a w mo d q)- 

To finish the simulation, we apply to the ancillas some (any) operator R that maps \aj) to 
hP0 9-1— J ). We then use a CNOT to copy the jth ancilla into the target qubit tj. We then undo 
all the previous computations to get rid of any conditional phase factors. As was remarked in 
Section|H3 the reverse computation uses W, which can be simulated exactly by evolving via G n 
for a positive time only with certain restrictions on the value of b (the value c is unimportant in 
that it only results in an overall phase factor). If q > 3, then it is easy to check that for real 
u > 0, exp(— iuG n ) oc / (restricted to the space of encoded vectors) if and only if both ua "^ +b ^ 
and nctn ( 2+b ) are both integers. The latter conditions hold iff b is rational. 

Finally, we note that we may be able to get by with less than the full use of E> and E on the 
inside of U and W in the circuit of Figure [TJ If the encoded state after applying U is not completely 
entangled, we need only decode the ancillae and those qubits that are entangle with the ancillas, as 
was done in the circuit of Figure |1J 
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5 Further Research 



We have assumed throughout that the coupling coefficients Jij of Q are all equal. Whether this 
assumption is realistic remains to be seen. It is certainly more likely in the short run that in feasible 
laboratory setups, the Jj.j will not be equal, but can still satisfy certain symmetries. For example, 
if n identical spin- 1/2 particles are arranged in a circular ring, we would expect the Hamiltonian 
to be cyclically symmetric, i.e., Jj j to depend only on (i — j) mod n. For another example, if the 
particles are arranged on points in a two- or three-dimensional regular lattice, we would expect 
translational symmetry of the Jij- 4 Computing parity in these more realistic situations would be 
very useful and deserving of further investigation. 

Heisenberg interactions also figure prominently in recent proposals for fault-tolerant quantum 
computation in decoherence-free subspaces (see, for example, |11[ IT] and references cited therein). 
The use of these interactions for this purpose does not appear consistent with our use here, yet it 
would be helpful to integrate these two approaches, perhaps by encoding logical qubits in a DFS. 
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